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Abstract The structure of the cosmic ray electron halo of 
a starburst galaxy depends strongly on the nature of galac- 
tic wind and the configuration of the magnetic field. We 
have investigated these dependencies by solving numerically 
the propagation of electrons originating in starburst galax- 
ies, most likely in supernova remnants. The calculations are 
made for several models for the galactic winds and for the 
configuration of the magnetic fields for comparison with ob- 
servations. Our simulation of a quasi-radio halo reproduces 
both the extended structure of ~ 9 kpc and the subtle hol- 
low structure near the polar region of the radio halo that are 
observed in the starburst galaxy NGC253. These findings 
suggest the existence of strong galactic wind in NGC 253. 
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starburst galaxies and the resultant radio halo to study the 
existence of the galactic winds and the configuration of the 
magnetic field. The calculation of the propagation of elec- 
trons is made by solving a coupled set of stochastic differ- 
ential equations (SDE) which is equivalent to the so-called 
diffusion convection partial differential equation and which 
has been successfully applied to the study of the solar mod- 
ulation phenomena of the galactic cosmic rays in the helio- 
sphere. In this paper we shall use spherically symmetric and 
axially symmetric models for the galactic winds and for the 
configuration of the magnetic field 1 1 2|. We present details 
of the structure of electron halo and also the resultant struc- 
ture of radio halo of starburst galaxies, with particular atten- 
tion to NGC 253. 



2 Numerical models 



1 Introduction 

Cosmic ray electrons accelerated in galaxies propagate out- 
ward by diffusion and convection in galactic winds and even- 
tually form an electron halo, the size of which depends on 
physical conditions such as the ambient magnetic field strength 
and strength of the galactic wind. The profile of nonthermal 
radio emission from the relativistic electrons reveals infor- 
mation on the nature of electron itself as well as informa- 
tion about the galactic wind from the starburst region and 
the configuration of the magnetic field. We have numerically 
investigated the structure of the cosmic ray electron halo of 
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The SDE equivalent to the diffusion convection partial dif- 
ferential equation is written as 



dX = udt + v/2KdW(t) , 

dP = -§P (V ■ u) dt - dP sync - dP IC , 



(1) 



where X and P are the position and the momentum of the 
particle, u is the galactic wind velocity, P syn c and Pic in- 
dicate the synchrotron and the inverse Compton momentum 
loss, k is the diffusion coefficient, and dW(t) is a Wiener 
process given by the Gaussian distribution. 
We adopted 
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where kb is the Bohm diffusion coefficient, B is the mag- 
netic field strength, ot is the Thomson cross section, F is 
the Lorentz factor of electrons, and U p h is the energy den- 
sity of the cosmic microwave background (CMB). Here we 
assume only CMB photons as the target photons for the in- 
verse Compton process. In our simulation, particles start at a 
fixed final point and run backwards in time until they come 
to the galactic disk boundary which has a 7 kpc radius and 
0.5 kpc thickness. The momentum spectrum /x(p) at arbi- 
trary position X is written as a convolution of the spectrum 
/x (po) at me galactic disk boundary with the normalized 
transition probability F(po, Xo|p, X) obtained by our SDE 
method as 

/x(p) = J fxo(po)F(po,X- \p,X)dpo . (5) 



Here the spectrum at the boundary is /x (Po) oc (m 2 c A + 
Po 2 c 2 )~ 16 /po which is assumed to be uniform at the galac- 
tic disk boundary. 

The simple analytical model for the spherically symmet- 
ric and the axially symmetric galactic wind flow from a star- 
burst galaxy are given by Chevalier and Clegg 1 1 1 and Zi- 
rakashvili and Voelk 1 2 1 respectively. We adopt their models 
for our simulations. 

For the spherically symmetric model of Chevalier and 
Clegg, the analytical solutions in the galactic halo are given 
as 



7 - 1 + 2/M 2 
1 + 7 



2\ (7+l)/[2(7-l)] 



R 



= ^ , (6) 



where u r and ug are the radial and latitudinal velocity com- 
ponent in cylindrical coordinates (r, 9, z), p is the density, p 
is the pressure, is the asymptotic velocity of the galac- 
tic wind, 7 is the adiabatic index, R g is the radius of the 
wind base, and p g is the p at r — R g . We adopted 7 = 5/3, 
R g = 300 pc and Uoa = 900 km/sec. 

We used the following expression |2| for the magnetic 
field which has a similar r-dependence as p(r, 9), 



B = B (sin6 



(13) 



where B is the magnetic field strength at the wind base and 
7 m is the adiabatic index of an isotropic random magnetic 
field. We assumed Bq = 50 pG and j m = 4/3. 



3 Results 



3.1 Structure of cosmic ray electron halo 
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where M is the Mach number, r is the radial coordinate, R 
is radius of the wind base, p is the density, P is the pressure, 
u is the galactic wind velocity, and 7 is the adiabatic index. 
We adopted R = 300 pc, 7 = 5/3, the mass loss rate M = 
3.9Af© yr _1 , and the energy production rate E = 1.9 x 
10 42 erg • s _1 . We also assume the magnetic field is frozen 

in the galactic wind flow, namely B{r) — (p(r)/ p(R)) 2 ^ 3 Bq, 
where Bo is the B at r = R. Bq is assumed to be 50 pG. 

For the axially symmetric model of Zirakashvili and Voelk, 
the analytical solutions in the galactic halo are given as 

U r = UooCOS^ , (9) 

ue = - (l - 7 _1 ) UooSin6» , (10) 

p M )=p s sin-^-j , (11) 

2-1 

p(r,9)= p g sm^9[^-) ^ \ (12) 



Fig. 1 (a) Side view of 1 GeV electron intensities integrated along the 
line of sight in the spherically symmetric model. The black rectangle 
indicates the galactic disk boundary. The peak intensity is normalized 
to 1.0. Contours indicate 1(T 6 , 10~ 5 , 10~ 4 , 10~ 3 , 0.01, 0.1, 0.5, and 
0.8.; (b) Same as (a) in the axially symmetric model. Contours indicate 
0.1,0.2, 0.3,..., 0.8, and 0.9. 

Fig-fflshows side views of the 1 GeV electron intensities 
integrated along the line of sight in the spherically symmet- 
ric and the axially symmetric models. The filled rectangle 
at the center indicates the side view of the galactic disk. 
We can see in Fig. [fl there is a big difference in the struc- 
ture and the size of the resultant electron halo of the two 
models for symmetry. One reason for this difference comes 
from the difference in the distribution of the source posi- 
tion of electrons on the galactic disk boundary. In the axi- 
ally symmetric model electrons may come from almost any 
point on the galactic disk boundary, however, in spherical 
symmetric model almost all electrons should originate near 
the galactic center where the magnetic field intensity is high 
and electrons suffer severe energy loss. The spherical sym- 
metry model clearly does not reproduce the observed radio 
halo of NGC253 |3 1. We discuss below only the results for 
axially symmetric models, because in many actual cases the 
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Fig. 2 Side view of electron intensities integrated along the line of 
sight in the axially symmetric model for 10 GeV electrons. The black 
rectangle indicates the galactic disk boundary. The peak intensity is 
normalized to 1.0. Contours indicate 0.1, 0.2, 0.3,..., 0.8, and 0.9. 
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synchrotron cooling time. For the axially symmetric model 
the TeV electron halo is not formed, because the synchrotron 
cooling time of TeV electrons is ~ 5 x 10 3 years, and the dif- 
fusion length is only about 1 pc. This is consistent with the 
H.E.S.S. result of TeV 7 -ray observation of NGC 253 H. 
The synchrotron energy loss is proportional to the square of 
magnetic field strength, therefore the weaker the magnetic 
field strength, the larger the size of the halo. Fig. [3] shows a 
side view of the 1 GeV electron intensities integrated along 
the line of sight in the axially symmetric model for a Bq of 
5 /itG. The size of the halo shown in Fig. [3] is much larger 
than that shown in Fig^> as expected. 

A slight hollow structure near the polar region in elec- 
tron halo is visible in Fig.^t and Fig.|3]for 1 GeV electrons 
and in Fig. [2] for 10 GeV electrons. This structure comes 
from the r dependence of energy loss in cylindrical coordi- 
nates (r, 9, z). Fig. shows the mean total energy loss of 
1 GeV electrons at z = 5kpc as a function of r. The en- 
ergy loss peaks near the polar region, as the energy loss is 
dominated by the adiabatic loss, as shown in Fig.0J>, by the 
convection due to the galactic wind and as the wind velocity 
is higher in polar region as shown by Eq[9] Fig. |4j> shows 
the r dependence of the fraction of the energy loss shown in 
Fig. 0Ji for the three processes, adiabatic, synchrotron, and 
inverse Compton. The slight hollow structure for 10 GeV 
electrons seen in Fig. |2] comes from the dominance of syn- 
chrotron loss over the other two energy loss mechanisms. 
(Due to page limitations we cannot present a similar figure 
to Fig.|4]in this paper.) The rate of synchrotron energy loss is 
proportional to the square both of energy and magnetic field 
strength. The magnetic field strength is higher in the polar 
region as shown by Ea ll3l Electrons which arrive at the po- 
lar region originate from near the galactic center where they 
suffer severe energy loss at the start of their journey. 



Fig. 3 Side view of 1 GeV electron intensities integrated along the line 
of sight in the axially symmetric model for a Bo of 5 /iG. The black 
rectangle indicates the galactic disk boundary. The peak intensity is 
normalized to 1.0. Contours indicate 0.1, 0.2, 0.3,..., 0.8, and 0.9. 



starburst region has an ellipsoidal shape or a disk geometry 

Propagation of electrons is governed by the diffusion 
process and convection in the expanding galactic wind and 
the associated energy loss processes. The diffusion process 
depends on the electron energy, however the convection does 
not. The adiabatic loss is proportional to the electron en- 
ergy, and the other energy loss processes are proportional 
to the square of electron energy. The synchrotron energy 
loss is larger than the inverse Compton energy loss, because 
B 2 /8ir ^> Uph- Accordingly the resultant structure depends 
on the energy of the electrons as we will see below. Fig. |2 
shows a side view of the electron intensities integrated along 
the line of sight in the axially symmetric model for 10 GeV 
electrons. The contours in Fig. |2| were drawn with the same 
steps in Fig.QJb for easier comparison. The higher the elec- 
tron energy, the smaller the size of the halo, due to the shorter 
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Fig. 4 (a) Mean energy loss of electrons arriving at z — 5kpc with 
an energy of 1 GeV as a function of r. (b) The fraction of energy loss 
for three processes corresponding to the mean energy loss shown in 
(a). The solid line indicates the adiabatic energy loss. The dotted line 
indicates the synchrotron radiation loss. The dot-dash line indicates the 
inverse Compton radiation loss. 



The differential electron energy spectrum varies with po- 
sition due to the variation of the total energy loss. Fig. [5] 
shows the simulated electron energy spectra at several se- 
lected positions. The modulated spectrum shifts from the 
spectrum at the galactic disk boundary to the low energy re- 
gion without changing shape while the energy loss is dom- 
inated by the adiabatic loss, namely d\ogE oc const, where 
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Fig. 5 Simulated electron energy spectra at several selected positions. 
The solid lines indicate the energy spectra in the electron halo. The 
dotted line indicates the energy spectrum in the galactic disk. 

E is the electron energy. But the spectrum steepens in the 
high energy region due to the synchrotron energy loss. The 
bending point of the spectrum appears as expected at the 
electron energy where the adiabatic loss rate equals to the 
synchrotron energy loss rate 0. 



3.2 Comparison with observation 

NGC 253 is a nearby edge-on starburst galaxy at a distance 
of 2.5 Mpc. The galaxy has two spiral arms, and a bar of 
projected length 7 kpc. It has an extended nonthermal radio 
halo of ^9 kpc. It has also extended X-ray halo, possibly 
related to the galactic wind. 

Fig-EDshows side view of the 1 GeV electron intensities 
multiplied by B 2 and integrated along the line of sight in the 
axially symmetric model. As the power of synchrotron radi- 
ation is proportional to electron intensity times the square of 
magnetic field strength, Fig.|6|may mimic the radio intensity 
profile. 

Our result for the quasi-radio intensity profile qualita- 
tively reproduces the observed large radio halo of ~ 9 kpc. 
The contours in Fig. [5| are drawn with steps proportional to 
the steps in Fig. 2 of Carilli et al. 1 3 1 for ease of compari- 
son of the results of our simulations with the observed radio 
halo. Looking closer at the observed halo, shown in Fig. 2 of 
Carilli et al. |3|, we can recognize a slight hollow near the 
polar region and prominence far above the disk. Our quasi- 
radio halo reproduce these two features as seen in Fig. |6| 
The hollow structure comes from adiabatic cooling of elec- 
trons by the higher velocity galactic wind near the polar re- 
gion as shown in Fig. QJ). Our results suggest the existence 
of a strong galactic wind in NGC 253. In our simulation, the 
cosmic ray sources are assumed to be distributed uniformly 
in the galactic disk. The nature of galactic winds and mag- 
netic fields in starburst galaxies will be revealed by further 
simulation experiments based on models which take into ac- 
count a more realistic distribution of cosmic ray sources in 
galaxies. 
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Fig. 6 Side view of electron (1 GeV) intensities times B 2 integrated 
along the line of sight in the axially symmetric model. The black rect- 
angle indicates the galactic disk boundary. The peak value is normal- 
ized to 1.0. Contours indicate 0.005, 0.01, 0.015, 0.02, 0.025, 0.0375, 
0.05, 0.0625, 0.075, 0.1, 0.125, 0.25, and 0.5. 



4 Conclusions 

We have examined numerically the structure of the cosmic 
ray electron halo of starburst galaxies. We confirmed that 
the formation of electron halo extending to ~ 9 kpc from 
galaxies is possible for the axially symmetric model of the 
galactic wind flow and the magnetic field configuration. Fur- 
thermore, a subtle hollow structure of radio halo near the po- 
lar region observed in NGC 253 was reproduced. This struc- 
ture may come from large energy loss near at the polar re- 
gion by either the adiabatic energy loss due to strong galactic 
wind or the synchrotron energy loss due to relatively strong 
magnetic field there. These findings suggest the existence of 
strong galactic wind in NGC 253. 
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